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I. INTRODUCTION 

Interacting particles rotating in a harmonic confine- 
ment has become an experimentally realizable quantum 
system which has attracted plenty of attention during the 
last decade. The experiments in semiconductor quan- 
tum dots in the presence of a magnetic field showed 
Hund's rul^ and properties related to the quantum Hall 
liquid^^. In increasing magnetic field the electrons, po- 
larized by the field, first form so-called maximum density 
droplef^, which is the finite electron number counterpart 
of the integer quantum Hall effect, and then are expected 
to form fractional quantum Hall liquid^ ^ ^ and eventually 
get localized to a Wigner crystal^. 

The weakly interacting boson and fermion gases in 
atomic traps are other systems where particles move 
in a nearly harmonic trap. The traps can be three- 
dimensional or quasi-two-dimensional. Rotating laser 
field has been used to put the atom cloud in rotational 
motion which has lead to observa tion of vortices in both 
boson and fermion system^^l^SI^. 

The experiments have inspired numerous theoretical 
studies on the properties of rotational quantum states 
of particles in a harmonic confinement. For quantum 
dots these studies have revealed several kinds of internal 
symmetry breakings of the system: Spin-density waveJ^^, 
broken symmet ry edgd^ , vortex format ioEp^Il^ and elec- 
tron localizatioiP^EZEHl^ In the case of boson and fermion 
condensates the theoretical research has concentrated on 
quite small angular momenta and has mainly been in- 
terested in vortex formation^^ althoug h local- 
ization by rotation has also been studiecP- ^ * ^^ * ^ I It has 
been observed that both the vortex formation^ and the 
localization^^ are qualitatively similar for both types of 
particles, and that the type of the interparticle interac- 
tion does not seem to play an important role. 

In this paper we will study in detail the limit of the 
particle localization to a 'Wigner molecule' where the 
many-particle spectrum can be determined by combining 
classical vibrational modes to the rigid rotation of the 
'molecule'. In the ca se of electrons this limit was studied 
in detail by MaksymF^ and later by uJ^. Here we will 



extend this work to bosonic systems and also study the 
effect of the range of the interparticle interaction. 

The paper is organized as follows. In Section II we in- 
troduce the models and discuss an exactly solvable model 
of harmonic interparticle interaction. In Section HI we 
present results, first for seven particles interacting with 
Coulomb interaction and then study the effect of the 
range of the interaction using as an example four particles 
interacting with a Gaussian interaction. The conclusions 
are given in Section IV. 



II. THEORETICAL MODELS 

A. Quantum mechanics in the lowest Landau level 

We assume a generic model of spinless fermions or 
bosons in a two-dimensional harmonic potential. The 
Hamiltonian is on the form 



H = 



2m 



N N N 

^V?+5:imc.o\^ + ^t;(|r,-r,|), (1) 



where N is the number of particles and m the electron 
mass and coq the oscillation frequency of the confining 
potential. The variable r = (x, y) is a two-dimensional 
position vector. The interparticle potential v is either 
the Coulomb interaction or a Gaussian repulsion: 



v{r) 



or v{r) 



(2) 



We solve the many-particle problem for a fixed angular 
momentum using the single particle basis of the lowest 
Landau level (LLL): 

where i is the single-particle angular momentum and 
a normalization factor. The restriction to the LLL is 
a good approximation when total angular momentum L 
is relatively large, L > N{N — ijMOZI^ corresponding to 
the filling factor of about u < 1/2 of the LLL. At this 
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point it is convenient to stress that, when we consider 
electrons, we do not exphcitly include magnetic field in 
our calculations. Assuming complete polarization, the 
only effect of the magnetic field has, is to increase the 
total angular momentum. 

We present our results in terms of the total angular 
momentum, which can be related to the filling fraction 
of the LLL hy u ^ N{N — 1)/2L. Angular momenta L = 
qN{N — l)/2 then correspond to the fractional quantum 
Hall liquid with v = 1/0, where q is the exponent of the 
Laughlin wave functioiP (Eq. (4) below). 

The restriction of the basis to the LLL gives some im- 
portant exact results about the properties of the quantum 
states. These can be seen by writing the Hamiltonian (1) 
in the occupation number formalism as 

H = Hq^XV = huJoN^^hcJoiclci^X^VijkiclcyiCk, 

£ ijkl 

(3) 

where cj creates a particle in the LLL. The states are 
written as \non1n2 • • • ) where the subsrcipt refers to the 
single particle angular momentum i. 

We notice immediately that for a fixed total angu- 
lar momentum L = n£ the diagonal term Hq of the 
Hamiltonian gives a constant huJo{L + A/"), i.e. it is inde- 
pendent of the configuration. The diagonalization of the 
Hamiltonian thus equals to diagonalization of the inter- 
action part V of the Hamiltonian, which gives directly 
the interaction energy AE = (^IV"!^). Moreover, this 
means that the structure of the quantum states are inde- 
pendent of coq and the coupling constant A, which have 
only the roles of giving the distance and energy scales. 

In the case of bosons the state with L = is the 
Bose condensate where uq = N and all the other states 
are empty, for example, for seven bosons the state is 
1 7000 • • • ) . In the case of fermions the lowest possible an- 
gular momentum is Lmdd = N{N — l)/2 and the state 
is the so-called maximum density droplet^ = 1 for 
i < N and this state for seven fermions is|111111100---). 
The only possible state at for seven bosons at L = 
1 is |61000---) and the corresponding state for seven 
fermions at L = Lmdd + 1 is |11111101000 • • • ). 

It follows that the interaction energy of the L = (or 
MDD) state equals to that of the L = 1 (or Lmdd + 1) 
state, for any two-body interaction. This is due to the 
fact that the L = 1 state is the center of mass excitation 
of the L = state. In the LLL the two lowest angu- 
lar momentum states for both bosons and fermions are 
thus completely fixed by the basis and, consequently, are 
independent of the interparticle interaction. 

In the case of the delta function interaction, the spin- 
less fermions do not feel the interaction due to the Pauli 
exclusion principle. The MDD can be written in the form 
of the Laughlin wave function (without normalization) as 

N 

= Yli^i - Zjye^'^''o\^^\'/^^, (4) 

i<j 



where z = x -\- iy, and the exponent q = I correspond- 
ing to the filling factor u = 1/q = 1. Since any wave 
function for higher angular momenta can be written as 
a symmetric polynomial times the above, the effect of 
the delta function interaction disappears. In the case of 
bosons the delta function interaction has an effect up to 
the angular momentum L = N{N — 1) when a symmetric 
wave function of the type of Eq. Q can be reached, with 
q = 2. 

It turns out that also for long range interactions, even 
for the Coulomb interaction, the Laughlin type wave 
function becomes rather accurate for g = 3 and 5 for 
fermions and for q = 2 and 4 for bosons, while at even 
higher values of angular momenta the particles localize 
to Wigner molecules which rigidly rotate in the harmonic 
well. 

While the Laughlin state gives good estimates for the 
wave function only for particular angular momenta, the 
composite fermion picture of Jain^^ gives estimates 
for all angular momenta. It is important to note that 
both Laughlin and Jain wave functions are derived from 
very general arguments and are completely independent 
of the interparticle interactions. 

The repulsive harmonic interaction 

^(|r,-r,.|) = -A|r,-r,f (5) 

is exactly solvabl^23^ESI^ xhe resulting interaction en- 
ergy of the lowest state as a function of the angular mo- 
mentum is for boson^^ 

AE{L) = -^N{N - 1) - ^NLO{L - 2) (6) 

where 6 is the step function. Due to the symmetry of the 
interaction the degeneracy of each state is large. Note 
that the interaction energy in this model decreases lin- 
early with L. 

B. Classical particles 

Classical electrons in a 2D harmonic potential ar- 
range in consecutive circles i.e. they form Wigner 
molecule d^^ l ^^l The molecule can have angular momen- 
tum due to the center of mass rotation or due to the 
rigid rotation of the molecule. In addition the molecule 
can have internal vibrational motion. The internal vibra- 
tional modes can then be found by solving the dynamical 
matrix. We are interested in the excitations of a rotat- 
ing system and therefore we have to solve the internal 
vibrations in a rotating frame of referen ce. For quan- 
tum dots this was first done by MaksymP^ who used 
the Eckart rotating frame. We use a more straigthfor- 
ward method where the Newton equations of motion are 
determined in a standard rotating frame, taking into ac- 
count the Coriolis force^^. We fix the rotation frequency 
cOr which determines the total angular momentum L of 
the rotating molecule. 
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The equilibrium positions of the particles depend nat- 
urally on the angular frequency and total angular mo- 
mentum L = IcOr^ where / is the moment of inertia 
I = m^rf . The classical energy 



1 



N 



i<j 



v{Ti 



21 



(7) 



is minimized in order to solve the equilibrium configu- 
ration. The interparticle interaction v{\ri — rj\) can be 
either Coulombic or Gaussian just like when solving the 
energy spectra of corresponding quantum mechanical sys- 
tems. 

In this article we report results for up to seven bosons 
and fermions. The eigenfrequencies can be solved ana- 
lytically (or numerically) from the equations of motion 
in the rotating frame: The Newton equations are lin- 
earized around around the equilibrium positions of the 
particles and the vibrational modes are obtained with 
matrix algebra. When the classical vibrational frequen- 
cies are solved we quantize the energies with canonical 
quantization^^ and get an estimate for the total quan- 
tum mechanical energy, which in general can be written 
as 



Eci(L) = ES(L) 



^ huJk{L){nk + ^) + hu;o{ncM + 1), 

k 

(8) 

where uJk{L) is the vibrational frequency determined in 
the rotating frame with angular momentum L, Uk = 
0, 1, 2, • • • , and the last term corresponds to the cen- 
ter of mass excitations. The symmetry or antisymmetry 
requirement of the total quantum state dictates which 
vibrational states are allowed for a given L. For this rea- 
son we have to determine the symmetry properties of the 
vibrational modes and use group theor}^'^^Eni3II3II to re- 
solve the allowed states. Using the notations of MaksynP 
we have the following requirements for fermions 



L + ^ riiki = I 



mod m for m = odd 
m/2 mod m for m = even. 



The corresponding requirement for both even and odd 
number of bosons is 



L + ^ riiki — mod m. 



The symbols refer to the number of excitations of 
mode uJi and the number ki to the symmetry property 
of the same mode (rotation by 27r/c/m is presented by 
exp(i27r/c/m))^. 

The classical system with a repulsive harmonic interac- 
tion can also be solved analytically. For the lowest energy 
state as a function of the angular momentum we get 



E{L) = L 



Quantizing this and taking into account the zero point 
energy, the interaction energy will be 



AE{L) =NuJoJl 



2{N -1)X 



2NX 



(10) 



2NX. 



(9) 



which for small A gives the result of Eq. ([6|. Notice 
that for cOq < 2NX the harmonic potential is not strong 
enough for confining the particles. 

We should also notice that in the case of the delta 
function interaction the classical system does not have 
any vibrational modes, but the particles behave like non- 
interacting particles (the collision probability is zero). 



III. RESULTS 
A. Seven particles with Coulomb interaction 

We have considered four and seven bosons and 
fermions with different interparticle interactions. We 
remind the reader that in the LLL the results are in- 
dependent of cjo, apart from the energy scale. It is 
then convenient to present the results in atomic units 
{m = h = e = ao = 1). We use ujq = 1 for the Gaussian 
interaction and cjq = 1/2 for the Coulomb interaction. 

The classical model has been studied for electrons 
in quantum dots by Maksym et al^^^, Matulis and 
Anisimovas^^ and by us^^. We will first demonstrate 
that the same model works also for bosons interacting 
with the Coulomb interaction. Since the earlier compu- 
tations have been restricted to six or less electrons we do 
this comparison for seven particles. 

The classical problem of seven charged particles in a 2D 
harmonic potential can be solved analytically by first de- 
termining the classical minimum energy ^ at the equi- 
librium in the rotating frame and then solving the dy- 
namical matrix. The lowest energy classical configura- 
tion is a hexagon of six particles surrounding the seventh 
particle at the origin. This geometry naturally has a six- 
fold rotation axis. The other stable classical geometry is 
a seven-fold ring. However, its energy is so much higher 
that it does not contribute to the low-energy rotation- 
vibration spectrum. 

The classical vibrational frequencies (except the center 
of mass mode) for 7 particles are shown in Fig. [l] The 
eigenfrequencies are shown as a function of the angular 
velocity of the rotation cjr- The angular momentum of 
the system increases rapidly and diverges as uor 
(0.5 in this case). Some values of the angular momenta 
are shown in the figure. We notice that all frequencies 
approach to either or 2iJo. O ne of the in ternal vibration 
energy is on the form c^bm = V^Sc^o^To^, and it is shown 
in the figure as a dotted line. This is a w ell known result 
for electrons forming; one single rin 
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FIG. 1: Upper panel: Classical vibrational frequencies as 
a function of the angular velocity of the rotating Wigner 
molecule with seven charged particles. The breathing mode is 
shown as a dotted line. The dashed and solid lines show other 
modes for seven electrons (in atomic units with uq = 0.5). 
Lower panel: Modes of the low energy vibrations (see text). 
Mode is similar to the mode cjs (and degenerate at Ur = 0) 
but the motion of the particles is opposite to that of mode 
(jjs. In mode 004 the particle at the center does not move. 
Mode ujr is similar to mode 004 (degenerate at Ur = 0) but 
the particles move opposite direction. Note that in mode coq 
the orbits of the particles are the same as in modes 004 and 
007 but the particles are at different phases. 



the breathing mode. The most interesting modes are the 
ones that approach 0. We have labeled these low-energy 
modes as ujs (lowest mode, plotted as red dashed line, 
ujs ^ 0.4 when Ur = 0), uo^ (dashed line, is degenerate 
with COS when cOr = 0), coq (solid line), 004 (red dashed 
line, 004 ^ 0.55 when cOr = and as C07 (the mode that is 
degenerate with uj^ when cj^ = and uj^ > UJ4 elsewhere) . 
We denote the center of mass vibrations as uoi and U02. 
The symmetry properties (i.e k values) are determined 
as iiP and they are for the low energy vibrational modes 
of seven electrons as follows: /ci = 5, /C3 = 1, /C4 = 4, 
/c5 = 5, /ce = 3 and kj = 2. We have visualized the 
electron motion in some of modes in the Fig. [T] 

Once we have the energies and the symmetry proper- 
ties of the classical vibrations solved, we need to quantize 
them to get the rotational spectrum of Eq. (8). De- 
termination of the symmetry properties means that for 
fermions we demand that the quantized state is antisym- 
metric and for bosons symmetric. 

Figure [2] demonstrates how the quantum mechanical 
spectrum approaches to that obtained by quantizing the 
classical vibrations. The figure shows the lowest (so- 
called yrast state) and second lowest energy for each an- 
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FIG. 2: Dependence of the ground state energy (solid lines) 
and the first excites state energy (dashed lines) on the total 
angular momentum in the case of seven particles interacting 
with the Coulomb interaction. The black lines (red lines) 
show the energy difference between the result of the exact di- 
agonalization (model Hamiltonian (8)) and the energy of the 
rigid rotation of the classical Wigner molecule. The upper 
panel is for bosons the lower panel for fermions. The inset 
shows the energy convergence as a function of the cut-off an- 
gular momentum of the LLL single particle basis in the case 
of fermions with L — 123. 



gular momentum. The energies are shown as the energy 
difference from the classically determined rotational en- 
ergy, Eq. (8) with Uk = ncM = 0. 

The results show the characteristic oscillation as a 
function of the angular momentum. The oscillation can 
be understood as a rigid rotation of the 'molecule' of lo- 
calized particles. The symmetry properties allow only 
every sixth angular momentum due to the six-fold sym- 
metry of the hexagon. For bosons the allowed angular 
momenta are L = 6n and for fermions L = 6n + 3, where 
n is an integer. For angular momenta in between, the 
purely rotational state has to be accompanied by a vi- 
brational state, and consequently the energy is higher. 

Figure [2] shows that already from angular momentum 
of about L = 40, corresponding to filling factor u = 1/2, 
the model Hamiltonian (8) describes qualitatively well 
the oscillations of the energy as a function of the angular 
momentum, when compared to the results of diagonaliza- 
tion of the Hamiltonian (3). The agreement gets better 
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with increasing angular momentum. Naturally, the spec- 
tra for fermions and bosons become then identical, apart 
of a phase shift dictated by the symmetry requirements. 

The reason that at high angular momenta the absolute 
energies of the exact diagonalization remains higher than 
those of the model Hamiltonian is due to the anharmonic- 
ity of the classical vibration potential: The zero point en- 
ergy contribution coming from the high energy modes of 
Fig. 1 causes the energy shift (i.e. the zero-point energy 
is underestimated due to the harmonic approximation). 
The restriction of the basis within the lowest Landau level 
or the cut-off used for large angular momenta do not play 
any marked role. The inset in Fig. 2 shows the depen- 
dence of the energy on the cut-off single particle angular 
momentum for L = 123. Note the different energy scale 
and that the total energy indeed depends very little on 
the cut-off angular momentum. 

Figure |3] shows in more detail part of the low energy 
spectrum for seven electrons calculated by two different 
ways, quantum mechanically and classically. We remind 
that in both cases the values are for interaction energy 
i.e = E — {N + L)hx)Q^ where E is either the solu- 
tion from the diagonalization of the quantum mechanical 
Hamiltonian (3) or the classical energy ([8|. In Fig. 3 
the classical spectrum is shifted up by a constant energy 
(0.014 a.u.) to make the comparison of the structure of 
the spectra easier. We can see that the two spectra agree 
in all details. Note also that the difference in the interac- 
tion energy given by the two models, 0.014/2.60 ^ 0.5%, 
is very small indeed. 



spectrum is obtained for fermions and bosons and the ap- 
proach to the classical limit is similar in both cases. 

Finally, we would like to note that in the case of seven 
particles the classical geometry of a seven-fold ring is also 
stable. However, its energy is so much higher that it does 
not have any contribution to the low energy spectrum at 
any angular momenta. The situation is different to that 
of six particles, where the ground state of a five-fold ring 
(with a particle at the center) and the six-fold ring are 
very clos e in e nergy and both contribute to the rotational 
spectrunP^. Nevertheless, at higher energies the seven- 
fold ring also exists. As an example. Fig. [4] shows the 
pair correlation function of the 17*^ state at L = 133. 
Clearly, it shows seven particles in a ring with an empty 
center. 



1 \ \ \ \ \ r 



_l 1 1 


1 1 i_ 








Jyj- 


~i 1 1 


{oJj - 

1 1 r 



FIG. 4: Pair correlation functions for selected states of 
seven electrons. Left panel shows the purely rotational state 
(ground state for L — 117). The center panel shows a vibra- 
tional state (fourth excited state for L — 133), which consists 
of vibrational modes cjs and cje- The right panel shows a 
seven-fold ring as a high excited state (17th state for L — 133). 




FIG. 3: Low energy spectrum of seven electrons at high an- 
gular momenta. Black dots show the exact diagonalization 
results and open squares the results of the model Hamilto- 
nian (8), shifted upwards with a constant 0.014. is the 
interaction energy in atomic units (a;o = 0.5). The ground 
states for L — 129 and 135 are purely rotational states. The 
ground states from L = 130 to L = 134 have vibrational 
modes c^s, (^6, ^7, ^3, respectively. 



The oscillatory behavior of the results presented in 
Figs. [2] and [3] is not surprising since similar results were 
reported by Maksym^ for three, four, five and six elec- 
trons in a quantum dot. Here we have shown that similar 



B. Gaussian interparticle interaction 

In order to study the effect of the range of the interpar- 
ticle interaction to the localization of particles we used 
Gaussian potential, Eq. (2). In this case we study a 
rather simple case of four particles. The ground state 
classical configuration of four particles is a square. Prus 
et aP^ have developed a model that considers electrons 
as Gaussian probability densities. They showed that this 
model yields also to the formation of a Wigner molecule. 
At this point we want to stress that our model is differ- 
ent from that because here the particles are point like but 
the interaction between the particles is Gaussian. Solving 
the eigenmodes is a bit more tricky than it was with the 
Coulomb interaction because the form of the interparti- 
cle interaction leads into transcendent equations at some 
point. As in the case of Coulombic interaction, knowing 
the modes is important, because the symmetry proper- 
ties of the mode define all the allowed combinations of 
angular momenta and vibrational modes. In the three 
particle situation this point is reached at the end of the 
calculations so that it is easy to keep accounts for which 
energy value belongs to which mode. The four parti- 
cle problem, however, is somewhat more troublesome be- 
cause the stage of transcendent equations comes before 
the diagonalization of the dynamical matrix. This leads 
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to difficulties of keeping track of eigenmodes and eigener- 
gies for different angular momenta. We solved this prob- 
lem by writing a program that finds out the right mode 
for a certain eigenenergy for every angular momentum 
separately. 

Due to the symmetry, the vibrational modes are iden- 
tical with the modes obtained in^^ with Coulombic in- 
terparticle interaction, and consequently their s ymm etry 
properties are known from the previous rep or t^^^. In 
the rotating frame the four particles have two low-energy 
modes. In the case of Gaussian interaction these two 
modes cross each other when the angular momentum is 
increased. The angular momentum where this happens 
depends on the width a of the Gaussian potential. This 
as such is interesting, and if the classical model describes 
the behavior of the quantum mechanical system correctly, 
the same phenomenon should also occur in the quantum 
mechanical spectrum. We will return to this detail later 
after we have presented the results for quantum mechan- 
ical systems and the corresponding classical systems. 




FIG. 5: Interaction energy for four fermions (uper panel) and 
four bosons (lower panel) as a function of the angular mo- 
mentum, calculated using different widths a of the Gaussian 
interaction. Black dots connected with solid lines show the 
results of the exact diagonalization in the LLL and the open 
circles connected with dashed lines the results of the model 
Hamiltonian (8). Red curves, a — 3.0; blue curves, a — 2.0; 
black curves; a = 0.5. The solid black and red lines (without 
points) show the exact and model Hamiltonian results for a 
harmonic interparticle interaction. 

The quantum mechanical and classical energies for 
Gaussian repulsion in the case of four fermions and 



bosons are represented in Fig. |5] for three different val- 
ues of cr: 0.5, 2 and 3. The results show that the 
behavior of the quantum mechanical systems approach 
to the behavior of corresponding classical systems when 
the Gaussian interparticle potential is broadened from a 
nearly contact potential into a smoother potential. When 
(7 = 3, the quantum mechanical spectrum is nearly iden- 
tical with the corresponding classical spectrum. When a 
is decreased the classical and quantum mechanical results 
start to depart, although the similar structure survives 
until very small values of <j. 

Figure |5] also shows the results for a delta- function in- 
teraction (cr = 0) and for a repulsive harmonic interac- 
tion {y{r) = l/(7rcr^)(l-r^/cr^), with cr = 3). In the case 
of the delta function interaction the interaction energy of 
the classical system is zero while the quantum mechanical 
particles feel each other up to Laughlin states where the 
contact interaction disappears, i.e. up to L = N{N—l)/2 
for fermions and L = N{N — 1) for bosons. Beyond these 
angular momenta, the quantum mechanical and classical 
energies trivially agree. At this point it is interesting 
to note that in quasi-one-dimensional quantum rings of 
fermions, where the spectrum is dominated by particle lo- 
calization at all angular momenta, even a delta function 
interparticle interaction leads to vibrational states^^. 

In the case of the repulsive harmonic interaction, Eq. 
(5), the classical result approaches fast the quantum me- 
chanical result when A increases, as seen from Eqs. (6) 
and (10). Figure [s] show results for the value A cor- 
responding to the leading term of Gaussian interaction 
with a = 3.0. In this case the classical and quantum 
mechanical results follow the same line shown in the fig- 
ures. Since the classical approach gives correct results in 
the delta function limit and for the wide Gaussian limit 
it is not surprising that it describes the quantum me- 
chanical spectra nearly exactly for all values of a when 
L > N{N -1). 

We have also made similar systematic study for three 
particles and limited computations for seven particles us- 
ing the Gaussian interactions. The results are qualita- 
tively similar to those presented here for four particles. 

The classical model seems to describe very well the 
quantum mechanical system especially at higher values 
of cr. In figure [6] we present a more detailed study of the 
effect of cr to the energy spectra. On the upper panel we 
present the ratio of classical and the quantum mechan- 
ical ground state energies as a function of cr for a few 
fixed angular momenta. On the other hand the results 
of the lower panel show how the radius of the classical 
equilibrium polygon (i.e a square) depends on the value 
of cr. The radii are normalized to one for every angular 
momenta when a ^ in order to draw the different radii 
on the same picture. The behavior of the radii means 
that for a fixed angular momentum the equilibrium con- 
figuration (i.e a square) ffist expands as a function of cr 
until it reaches its maximum value, and after that starts 
to shrink. The results on the two figures show a signifi- 
cantly similar behavior of the two variables as a function 
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FIG. 6: Upper panel: Ratio of the ground state energy of 
the classical model Hamiltonian and the exact diagonaliza- 
tion energy Ed/Ec^ as a function of the width of the Gaussian 
interaction a. Lower panel: Radius r of the classical equilib- 
rium polygon (a square) as a function of cr. ro is the radius 
for a — 0). Dashed lines are for angular momentum L — 6, 
solid lines for L = 10 and the thick solid lines for L = 18. 



of <j. The explanation of such similarity is that classical 
particles do not feel each others at all if the Gaussian 
potential is steep. When the potential broadens the po- 
tentials of different particles start to overlap. Quantum 
mechanical particles on the other hand feel each others 
even if the interparticle potential is a point contact (i.e 

cr ^ 0). 

We will now return to the crossing of the energies of the 
two lowest vibrational modes, when the angular momen- 
tum increases. The symmetry dictates that one of these 
modes is the lowest energy state at angular momenta 
L = 4n in the case of fermions (L = 4n + 2 for bosons) 
and the other mode at angular momenta L = 4n + 1 
for fermions (L = 4n + 3 for bosons). In order to get 
an estimate of the vibrational energy from the exact di- 
agonalization we have to subtract from the total inter- 
action energy the energy of the rigid rotation. This we 
have estimated by fitting a smooth function of angular 
momentum through the energies of the purely rotational 
states L = 6, 10, 14, • • • for fermions. The resulting vi- 
brational energies for the Gaussian width a = 3.0 are 
plotted in Fig. [?) together with the classical vibrational 
energies, as a function of the angular momentum. The 
results show that both classical and quantum mechanical 
vibrational energies cross each other at about the same 



0.006 



/ 


\^ 




\ ^ ^ 


\/ 


^^^^ \ 


/\ / 




/ \/ 


^ ^ ^^^^ ^ 



10 15 20 25 30 35 40 45 50 

L 



FIG. 7: Dependence of the two low-energy vibrational modes 
on the total angular momentum in the case of four fermions 
interacting with a Gaussian interaction with a — 3.0. Dashed 
red lines show the classical results and solid black lines the 
quantum mechanical results from the exact diagonalization. 



angular momentum. This finding supports the previous 
notions that the properties of the quantum mechanical 
system of a few particles rotating in a harmonic confine- 
ment can be explained with finest details by a model of 
vibrating molecule of localized particles. 
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FIG. 8: Average radius of the ring of four bosons as a function 
of the angular momentum L, calculated with exact diagonal- 
ization in the LLL. The center line and points show the radius 
R and the lower and upper lines and points VK, where W 
is the variance of the radius. Solid lines show results for the 
Gaussian interaction with g — 3.0, squares for a — 0.05 and 
crosses for the Coulomb interaction. Dotted line is the result 
for classical particles with Gaussian interaction, cr = 3.0. 

Finally, we would like to quantify the degree of local- 
ization in the quantum mechanical many-particle system. 
The symmetry of the Hamiltonian dictates that the to- 
tal particle density is circularly symmetric. However, for 
small number of particles they localize in circles and con- 
sequently the radial density distribution gives direct in- 
formation of the degree of localization. In the case of four 
electrons we have determined the average radius R of the 
ring of electrons and its variance W, which describes the 
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FIG. 9: Pair correlation functions of four fermions interacting 
with the Gaussian interaction with a — 3.0. Upper row shows 
results for L = 10 (left) and L = 18 (right) and the lower row 
for L = 30 (left) and L — 42 (right). The contour plots are 
in the same scale to demonstrate the expansion due to the 
rotation. 



creased momentum of inertia. For small number of par- 
ticles this happens already at small angular momenta. 

As an example we have studied the case of seven par- 
ticles interacting with Coulomb interaction. The ground 
state geometry of a nonrotating system is a pentagonal 
bipyramid, i.e. a five-fold ring in the x-?/-plane with two 
atoms at positions ±cz. When this is put in rotation 
the five-fold ring will expand while the particles in the 
z axis will get closer to each other. However, already at 
angular momentum L = 7 the planar hexagon with one 
particle at the center will have lower energy and at angu- 
lar momentum L = S the energy difference between the 
3D and 2D structures is much larger than the vibrational 
energies of the 2D structure. 

Similar results are obtained for other small particle 
numbers in 3D harmonic oscillator. Consequently, we 
can safely state that the region where the classical model 
explains the spectrum for the 2D harmonic oscillator the 
result would be the same had we done the computations 
for the 3D confinement. 



IV. CONCLUSIONS 



width of the localized electron state: 

= ^ / dvrpir), W' = ^l dv{r - Rfp{r), (11) 

where p is the total electron density. It turns out that for 
the quantum mechanical system both the radius and its 
variance are nearly independent of the shape or strength 
of the interparticle interaction. The results for bosons 
are shown in Fig. [8]as a function of the angular momen- 
tum. The radius reaches fast the classical value and it is 
nearly independent of the interparticle interaction. Simi- 
larly the width is independent of the interaction and also 
nearly constant as a function of the angular momentum. 
The increased localization is mainly due to the increase 
of the radius with L which makes the relative width to 
decrease. 

The same result is seen qualitatively in Fig. [9] where 
the pair correlation function for four fermions with Gaus- 
sian interaction (<j = 3.0) is shown for four different val- 
ues of the angular momentum. The pair correlation func- 
tion demonstrates that the angular localization increases 
with the same rate as the radial localization. 



C. Rotation in a 3D harmonic confinement 

Above we have considered rotational spectra of inter- 
acting particles in a 2D harmonic potential. In the 3D 
case the classical particles interacting with repulsive in- 
ter at ion tend to form spherical shells^^. However, when 
put into rotation, it is energetically favorable for the 3D 
structure to collapse into a 2D structure due to the in- 



We have studied the rotation induced localization 
of particles, in a two-dimensional circular harmonic 
trap. Earlier studies have already shown that a long- 
range Coulomb interaction causes formation of Wigner 
molecules in quantum dots and that the quantum me- 
chanical many-particle spectra can be determined by 
canonical quantization of the rigid rotation and vibra- 
tional modes of the molecule^ The present calcula- 
tions show that Coulomb interaction also causes simi- 
lar localization of bosonic particles. For seven particles 
(fermions or bosons) the low energy rotational spectrum 
is determined by vibrational modes of the molecule al- 
ready before the angular momentum corresponding to 
the filling factor u = 1/2 of the quantum Hall liquid is 
reached. 

For mimicking atoms in a harmonic traps we have used 
Gaussian interparticle interaction and studied the depen- 
dence of the localization on the range of the interaction. 
It was observed that both in the short range (delta func- 
tion) and long range (repulsive harmonic) limit the clas- 
sical model gives exactly the quantum mechanical ener- 
gies. Consequently, for any range of Gaussian interaction 
the classical model describes very accurately the quan- 
tum mechanical many-particle energy spectrum at large 
enough angular momenta. However, the delta function 
interaction does not support particle localization (parti- 
cles do not feel each other). 

We also estimated that in a three-dimensional har- 
monic confinement the rotating particles form very fast 
a two-dimensional configuration. For example, for seven 
particles with Coulomb interaction this happens already 
at L = 7, meaning that in the 3D confinement the low- 
energy rotational spectrum agrees with that calculated 
for a 2D system. 
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In conclusion, independent of the form of the repulsive 
interparticle interaction or the type of the particles, the 
rotational spectrum of the particles confined in a har- 
monic trap can be estimated by quantizing the classical 
vibrational modes determined for localized particles in a 
rotating frame. 
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